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Abstract 

We present results of the first fully-dynamical lattice QCD determination of hyperon-nucleon 
scattering. One s-wave phase shift was determined for nA scattering in both spin-channels at 
pion masses of 350, 490, and 590 MeV, and for nS~ scattering in both spin channels at pion 
masses of 490, and 590 MeV. The calculations were performed with domain-wall valence quarks 
on dynamical, staggered gauge configurations with a lattice spacing of 6 ~ 0.125 fm. 
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A. Introduction 



In high-density baryonic systems, the large value of the Fermi energy may make it energet- 
ically advantageous for some of the nucleons to become hyperons, with the increase in rest 
mass being more than compensated for by the decrease in Fermi energy. This is speculated to 
occur in the interior of neutron stars, but a quantitative understanding of this phenomenon 
depends on knowledge of the interactions among the hadrons in the medium. In contrast 
to nucleon-nucleon (NN) interactions, where the wealth of experimental data has allowed 
for the construction of high-precision potentials, the hyperon-nucleon (YN) interactions are 
only very-approximately known. Experimental information about the YN interaction comes 
mainly from the study of hypernuclei [2|, the analysis of associated A-kaon and E-kaon 
production in NN collisions near threshold jl, HJs], 0, 0, S], and hadronic atoms There 
are a total of 35 cross-sections measurements [lO|] of the processes Ap Ap, E~p — An, 
S+p — s> S+p, T,^p Tj^p and S^p S°n, and unsurprisingly, the extracted scattering 
parameters are highly model dependent. The theoretical study of YN interactions is hin- 
dered by the lack of experiment al g uidance. The "realistic" potentials developed by the 
Nijmegen 0, [ij and Jiilich [H, [U0 groups are just two examples of phenomenological 
models based on meson exchange. These are soft-core potentials with one-boson exchange 
models of the NN interaction. Since SU{2)) flavor symmetry is broken by the differences in 
the quark masses, the corresponding couplings are not completely determined by the NN 
interaction and are instead obtained by a fit to the available data. In Ref. [ll. 



12 , for exam- 



ple, six different models are constructed, each describing the available YN cross-section data 
equally well, but predicting different values for the phase shifts. The effective field theory 



approach [16|, |17|, Il8|, Il9|, [20 



is less developed and suffers from a large number of couplings 
that need to be fit to the data. 

In view of the large uncertainties in the YN scattering amplitudes and their importance 
for modeling neutron stars and the study of hypernuclei, a first-principles QCD calculation 
of YN scattering is highly desirable. The only way to achieve this is through numerical 
calculations using lattice QCD. In a previous paper [191], some of the present authors outlined 
a program to address this issue with a combination of lattice calculations and the use of 
effective field theories. This paper reports on the first results of the implementation of this 
program. In particular, we compute low-energy s-wave phase shifts for YN scattering in 
the ^5*0 channel and '^Sj_ 



^Di coupled-channels at particular energies, using Liischer's finite- 
This is a straightforward extension of work by some of the 
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volume method 

present authors on NN scattering from lattice QCD [2j]. The calculations were performed 
in the isospin limit on the coarse MILC lattices [2^ for pions with masses of ~ 290 MeV, 
~ 350 MeV, ~ 490 MeV and ~ 590 MeV. However, we do not attempt to extrapolate to the 
physical pion mass as it is likely that all but one of the data points lies outside the regime 
of applicability of the YN EFT's. 



B. Methodology and Details of the Calculation 



Lattice QCD calculations of the interactions among hadrons are notoriously difficult and 
require circumventing the Maiani- Testa theorem 2^, which states that one cannot compute 
Green's functions at infinite volume on the lattice and recover S-matrix elements except 
at kinematic thresholds. The s-wave scattering amplitude for two particles below inelastic 

which entails a measure- 



thresholds can be determined using Liischer's method [21|, |22|, [23 
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ment of one or more energy levels of the two-particle system in a finite volume. Our compu- 

laces 



tation uses the mixed-action lattice QCD scheme developed by LHPC [27|, |28[ which t 
domain-wall valence quarks from a smeared- source on Nf = 2 + 1 asq t ad- improved [2 
MILC configurations generated with rooted staggered sea quarks [25| that are hypercubic 
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smeared (HYP-smeared) [3l| 



In the generation of the MILC configurations, the 
strange-quark mass was fixed near its physical value, brris = 0.050, (where b = 0.125 fm is the 
lattice spacing) determined by the mass of hadrons containing strange quarks. The two light 
quarks in the configurations are degenerate (isospin-symmetric) , with masses bmi = 0.007, 



0.010, 0.020 and 0.030. As was shown by LHPC [27|,|28|, HYP-smearing allows for a signif- 
icant reduction in the residual chiral symmetry breaking at a moderate extent Ls = 16 of 
the extra dimension and domain- wall height M5 = 1.7. Using Dirichlet boundary conditions 
we reduced the original time extent of 64 down to 32. This allowed us to recycle propaga- 
tors computed for the nucleon structure function calculations performed by LHPC. For bare 
domain-wall fermion masses we used the tuned values that match the staggered Goldstone 
pion to few-percent precision 



27, 28 



can be found in Ref. 
software suite 36, 37 
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The parameters used in the propagator calculation 
All propagator calculations were performed using the Chroma 
. Eight propagators per configuration were computed at distinct source 
points on the lattice. 

We found that the cleanest method for extracting the energy- difference between the YN 
state, and the mass of an isolated nucleon and an isolated hyperon, was by forming the ratio 
of correlation functions 



G 



Y"Af(^) 



-AEyNt 



(1) 



where S denotes spin. At large times, this ratio depends exponentially upon the ground-state 
energy shift of the YN system due to interactions. The single nucleon correlator is 



CN{t) = N\0,0)) 



(2) 



and the single hyperon correlator has an analogous form. The YN correlator that projects 
onto the s-wave state in the continuum limit is 



Y.{Yr{t. ^)N^it, y)r;t(o, o)ivft(o, o)) 



(3) 



x,y 



where a,P,a,p are isospin-indices and i,j,k,l are Dirac-indices. The tensor X]^^lp has 
elements that produce the correct spin-isospin quantum numbers for a hyperon and nucleon 
in an s-wave. The summation over x (and y) corresponds to summing over all the spatial 
lattice sites, thereby projecting onto the momentum p = state of each particle seperately. 
The interpolating field for the proton is pi{t,x.) = eafecw"(t,x) (^u''^ {t,x.)C'y5d'^{t,x.)^ , and 
similarly for the neutron and hyperons. We have used an interpolating field n x S~ to 
determine the energy-eigenvalues of the s-wave strangeness = 1, isospin = | eigenstates in 
both spin channels, and an interpolating field n x A to determine the energy-eigenvalues of 
the s-wave strangeness = 1, isospin = | eigenstates in both spin channels. 

Once the energy shift due to the YN interactions has been computed, the real part of the 

To extract 



inverse scattering amplitude is determined via the Liischer formula [21|, |22|, [23 



pcot(5(p), where S{p) is the phase shift, the magnitude of the center-of-mass momentum. 
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is extracted from the energy shift, AEyn = yP"^ + My + yp"^ + Mj^ — My — M^, and 
inserted into: 

Pcot^(p) = ;^s(g) . (4) 
which is vahd below the inelastic threshold. The regulated three-dimensional sum is 

S(^) - - 4vrA , (5) 

j IJI / 

where the summation is over all triplets of integers j such that | j | < A and the limit A — oo 
is implicit. 



C. Results 



The effective mass plots of the ratio of correlation functions with identifiable plateaus ob- 
tained at rriTr ~ 350 MeV, m,r ~ 490 MeV and m^r ~ 590 MeV are shown in figs. [1], [2l 
|3] and HI Single and double exponential forms were fit to the correlation functions by 
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FIG. 1: Effective mass plots for nA in the ^S'o-channel at pion masses of m-,^ ^ 350 MeV (left 
panel), m-n ~ 490 MeV (center panel) and mjr ~ 590 MeV (right panel). The straight line and 
shaded region correspond to the extracted energy shift and associated uncertainty. The dashed 
lines correspond to the statistical and systematic errors added linearly. 



X^-minimization, from which the YN energy shifts were determined. The central values and 
statistical uncertainties were determined by the jackknife procedure over the ensemble of 
configurations, and are shown in Table [B 

The plateaus in the correlator ratios GyAr(t) persist for only a small number of time-slices. 
At small t there is the usual contamination from excited states whereas at larger t the signal- 
to-noise ratio degrades exponentially with t according to e~'^*^^"'"^'*'~^™'")*[38']. The Dirichlet 
boundary at t = 22 introduces a systematic uncertainty due to backward propagating states. 
However, in practice, the statistical noise becomes a limiting factor at far earlier time slices 
and the boundary at t = 22 is not an issue for this calculation. We obtained a non-zero 
energy shift larger than the statistical error in ten of the Gyn correlation functions, as shown 
in Table [H The phase shifts 5 (and — l/fccot5(fc)) were determined through the Liischer 
formula and their errors by the jackknife procedure (we do not give a value of 5 for a negative 
energy shift). The quantities in TableUthat are in physical units were obtained with a lattice 
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FIG. 2: Effective mass plots for nA in the ^5i-channel at pion masses of m,r ^ 350 MeV (left 
panel), ^ 490 MeV (center panel) and m.^ ~ 590 MeV (right panel). The straight line and 
shaded region correspond to the extracted energy shift and associated uncertainty. The dashed 
lines correspond to the statistical and systematic errors added linearly. 
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FIG. 3: Effective mass plots for nYr in the ^S'o-channel at pion masses of 771,7^ ~ 490 MeV (left 
panel) and txIt^ ~ 590 MeV (right panel). The straight line and shaded region correspond to the 
extracted energy shift and associated uncertainty. The dashed lines correspond to the statistical 
and systematic errors added linearly. 



spacing of 6 = 0.125 fm set by MILC, which is consistent with the determination from 
(6 = 0.1274 ±0.0007 ±0.0003 fm) based on the chiral expansion of the pion decay constant. 
We have not shown the results for channels in which there is no clear plateau in the effective 
mass plot. This is the case for all ~ 290 MeV {mi = 0.007) correlation functions where 
the rapid decrease of the signal-to-noise ratio caused by the small pion mass eliminated all 
plateaus. The systematic errors shown in Tabled are determined by varying the fitting range, 
and by comparing the results of fitting one and two exponentials to the ratio of correlation 
functions. 

It is not clear that we have been able to identify the ground states in all of the correlation 
functions, e.g. nS" in the ^S'o-channel at ~ 490 MeV, and nA in the ^S'o-channel at 
m^r ~ 490 MeV, as the statistics are not sufficient to determine whether the large-time 
behavior we observe is due to noise or due to the presence of any states with lower energy 
than those shown in Table [H Indeed, it would be very exciting if there were states with 
lower energy, as they would likely be bound states (based on naturalness arguments and the 
exact Liischer relation). This uncertainty in no way undermines our results; regardless of the 
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FIG. 4: Effective mass plots for riT," in the ■^S'l-channel at pion masses of rriT^ ~ 490 MeV (left 
panel) and rriTj ^ 590 MeV (right panel). The straight line and shaded region correspond to the 
extracted energy shift and associated uncertainty. The dashed lines correspond to the statistical 
and systematic errors added linearly. 



Channel 


(MeV) 


Range 


AE (MeV) 


|k| (MeV) 


6 (degrees) 


-(A; cot (5) -1 (fm) 


nA 


592 ± 1 ± 10 


8-12 


-9±8±20 






0.8 ± 1.4 ±0.4 




493 ± 1 ± 8 


6-9 


29.8 ±5.4 ±2.5 


197 ± 24 ± 4 


-32.3 ±8.1 ±2.8 


0.63 ±0.12 ±0.014 




354 ± 1 ± 6 


5-9 


56.8 ±6.0 ±5.5 


255 ± 22 ± 13 


-53.4 ± 8.5 ± 10.1 


1.04 ±0.24 ±0.15 


nA 


592 ± 1 ± 10 


8-13 


-13± 13±8 






3± 14±2 


's^ 


493 ± 1 ± 8 


7-11 


-4± 13± 14 






{—oo, oo) 




354 ± 1 ± 6 


5-10 


23± 17±4 


168 ± 62 ± 14 


-23± 18±4 


0.50 ±0.26 ±0.06 




592 ± 1 ± 10 


9-13 


-17± 11 ± 27 






(— oo, c») 




493 ± 1 ± 8 


5-9 


24.9 ±7.8 ±3.0 


179 ±28 ±11 


-27.2 ±9.0 ±3.8 


0.57 ±0.13 ±0.05 




592 ± 1 ± 10 


6-10 


38.5 ±8.8 ±5.0 


226 ± 26 ± 15 


-44.3 ±9.8 ±5.4 


0.85 ±0.20 ±0.10 




493 ± 1 ± 8 


6-10 


53 ± 14 ± 5 


261 ± 35 ± 13 


-58 ±15 ±5 


1.19 ±0.51 ±0.15 
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TABLE I: Summary of results from the Gyn correlation functions which exhibit a clear plateau 
in the effective energy plot. The first error is statistical and the second error is systematic. 

nature of the states shown in Table [H they are clearly YN states present in the continuum. 

In addition to the fitting systematics given in Table [U there are other systematic uncer- 
tainties in our calculations that we have not shown, as they are all expected to be small 
in comparison. The discretization errors due to the finite lattice spacing arising from the 
sea-action are of order O^asb"^) {as is the strong coupling constant) and those in the valence 
sector of order 0{h'^) due to the near-perfect chiral symmetry. The agreement between con- 
tinuum chiral perturbation theory and other results based on the same discretization scheme 
M H S, 13 strongly suggest that the discretization errors in the hyperon interactions 



are, at most, of the order of a few percent. This is much smaller than the statistical errors 
quoted in Table [H A similar point can be made regarding the use of different fermion actions 
in the valence and sea sectors. Further, the relation between energy levels and phase shifts, 
eq. (jlj), is valid only up to corrections that are exponentially small in the volume. The 
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corrections to the Liischer formula can be computed in chiral perturbation theory, as shown 
in the vrvr case in 41 1 and for two nucleons in j42|. These effects are particularly small in 
the NA system, as the long-range part of the interaction is dominated by two-pion exchange 
and one-kaon exchange, and not one-pion exchange. 
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FIG. 5: Comparison of the lowest-pion-mass lattice results in each channel with a recently develo ped 
YN EFT [2(J] (squares), and several potential models: Nijmegen n| (triangles) and Jiilich 15 1 
(diamonds). The dark error bars on the lattice data are statistical and the light error bars are 
statistical and systematic errors added in quadrature. 



D. Discussion 



We have presented results of the first fully-dynamical lattice QCD calculation of YN inter- 
actions. The scattering amplitudes for s-wave nA and nS~, in both the ^^o-channel and the 
^Si —^Di coupled-channels, have been determined at one value of momentum for pion masses 
of ~ 350 MeV, ~ 490 MeV and 590 MeV. Unfortunately, the lightest pion mass at which 
we have been able to extract a signal is at the upper limits of the regime of applicability 
of the effective field theories that have been constructed, thus precluding a chiral extrap- 
olation. However, this work does provide new rigorous theoretical constraints on effective 
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field theory, and potential model constructions of YN interactions. In fig. [5] we compare the 
lattice values of the phase shifts to recent EFT results [20| (squares), and to several potential 



models: Nijmegen [ll| (triangles) and Jiilich ISj. At face value these results appear quite 



discrepant, however one should keep in mind that extrapolation to the physical pion mass 
will seriously alter individual contributions to the YN interaction. 

While the measurements of the momenta and phase shifts are unambiguous, their physical 
interpretation is not entirely resolved. Each of the phase shifts at the lowest pion masses are 
negative. Assuming that the observed state is the ground state in the lattice volume, this 
implies that the interactions are all repulsive. The interaction in the ^5*1 —^Di coupled 
channels is strongly repulsive at m.,^ ~ 490 MeV, while the interaction in the ^^o-channel is 
only mildly repulsive. The opposite is found to be true for the riA systems at ~ 350 MeV, 
where the interaction in the ^^o-channel is found to be strongly repulsive, while that in 
the ^5*1 — ^-Di coupled channels is mildly repulsive. However, there may be channels for 
which there exist states of lower, negative energies, some of which may correspond to bound 
states in the continuum limit. If such states are present, then we would conclude that 
the interaction is attractive, and that the positive-shifted energy state we have identified 
corresponds to the first continuum level. Current statistics are sufficiently poor that nothing 
definitive can be said about the existence of such states. Therefore, we are continuing to 
accumulate statistics and experiment with signal optimization in order to resolve this issue. 

It is clear that a precise lattice QCD calculation of YN scattering will have dramatic 
impact upon the field of hypernuclear physics, and may have an equal impact on our ability 
to determine the evolution of neutron stars, simply due to the present absence of precise 
experimental data. We have performed the first of such calculations, albeit at unphysically 
large pion masses. The present work was limited entirely by the lack of computational 
resources. We hope that this limitation recedes in the future, and that lattice QCD can be 
developed as a reliable tool to calculate the interactions between baryons in experimentally 
inaccessible or difficult areas of strong interactions. 
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